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Abstract 

A one dimensional lattice model is formulated to study tapping dynamics and the long time 
steady distribution in granular media. The dynamics conserves the number of particles in the 
system, and density changes are associated to the creation and destruction of empty sites. The 
model is shown to be consistent with Edwards thermodynamics theory of powders. The relationship 
with lattice models in which the number of particles is not conserved is discussed. 
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I. INTRODUCTION 



In the last years, a great deal of effort is being made trying to understand the physical 
mechanisms leading to compaction in weakly vibrated granular systems, and the properties 
of the steady state eventually reached in the long time limit. This has been prompted and 
stimulated by the seminal papers of the Chicago group reporting experimental results of 
compaction [1|, 

as 

. Granular compaction consists in the increase of the density, starting 
from an initial low-density state, as a consequence of external excitations, usually vertical 
shakes or taps. Every tap is followed by a free relaxation, so that in the process the system 
goes through a series of blocked configurations. 

Starting from an "ergodic hypothesis" for powders, based on the extensive, global, char- 
acter of the dynamics induced by shaking, Edwards and coworkers have formulated a 
microscopic theory for the steady state of vibrated granular media that is similar to conven- 
tional statistical mechanics. Moreover, they assume that the steady state is fully character- 
ized by the volume of the system, which plays then a role analogous to that of the energy 
in usual thermal systems. This provides the "microcanonical" description. The associated 
"canonical" probability distribution is obtained by maximizing the statistical entropy under 
the condition that the average volume is given. Of course, the probability of a given con- 
figuration depends only of its volume. The parameter conjugated of the volume, similar to 
the thermal temperature, was named compactivity by the authors in Ref. ^J. 

Up to now, there is not a definite experimental test of the above thermodynamic theory 
of powders. The measure of the compactivity, or the entropy, of a granular system seems a 
rather difficult task not only in real experiments but also in realistic models, although some 
procedures have been proposed. They are based on the determination of the average volume 
and its fluctuations as a function of the control parameter of the shaking, e.g. the vibration 
intensity ^, 0|. From these two functions, the compactivity can be obtained, in principle, 
by integration, although this program is hard to carry out in practice due to the uncertainty 
of the measurements. Another alternative way, this one based on the generalization of the 
Einstein relation between diffusivity and mobility, has been recently discussed and analyzed 
in a system of inelastic hard spheres by means of molecular dynamics simulations [{J. 

On the other hand, the validity of Edwards theory has been studied in the context of 
several simple models, with different underlying physical mechanisms. It has been found 
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that the results for Tetris and spin-glass models are consistent with the theory. In 

these systems, the number of particles is fixed but most of the results are numerical, due to 
the complexity of the models used. One dimensional Ising models, with and without kinetic 
constrains, have also been considered jloL 14], because they are simple enough 

as to allow a detailed analytical study in many cases. For weak tapping, agreement with 
Edwards's theory was found again, although discrepancies show up in the limit of strong 
tapping. Quite interestingly, all the Ising-like models in Refs. |l0, III) 12, Q Q have 
been formulated as open systems. The number of particles does not remain constant, but it 
changes along the compaction process, as a consequence of adsorption-desorption events from 
a theoretical particle reservoir in contact with the system. Instead, it is the volume what is 
kept fixed, leading in this way to the variation of the density. Then, although it is true that 
the steady distribution of these models can be considered as a "grand canonical" ensemble 
generalization of the theory, it is also clear that it is not characterized by the compactivity 
(temperature) but by another parameter playing the role of the chemical potential. This 
difference is evidently relevant when trying to relate any of them with the characteristics 
of the vibration process, e.g. its intensity. Beyond that, the distinction might become 
conceptually crucial when dealing with granular mixtures and segregation phenomena. In 
that case, each of the different species is going to have its own analogous of the chemical 
potential parameter . Whether or not it is also needed to consider different compactivities 
(temperatures), as it as been suggested recently jlra j. is a different question. 

The aim of this paper is to present a closed, constant number of particles, one-dimensional 
lattice model for compaction. Again the model is simple enough as to be analytically 
tractable. During the tapping process, particles diffuse and also empty sites (holes) are 
created and destroyed in the system, according with well defined rules. The latter are chosen 
to mimic, in a crude way, what happens in real compaction experiments. More precisely, the 
model tries to represent a vertical section of a vibrated two-dimensional system. In a shake, 
the length of the section can increase because empty regions (holes) are created between 
particles. These regions can be used for the particles to diffuse. Afterwards, once the shake 
has ended, the system tries to compact due to the action of gravity. This is accomplished in 
the model by means of the elimination of holes. But, because of the geometrical constrains 
following from excluded volume effects of the hard particles in the neighbouring vertical 
sections, not all the holes can be destroyed in the free evolution. Only large enough empty 



3 



regions can be reduced. As a consequence of the combination of a tap and the next free 
relaxation, the length can increase in some regions of the system and decrease in others. 
The net balance determines the global behaviour of the system in the compaction process. 

The plan of the paper is as follows. In Sec. |H1 the model is formulated at the mesoscopic 
level of description by means of a master equation for the transition probability. This equa- 
tion is exactly solved for the steady distribution in Sec. IIII[ and the associated macroscopic 
description is discussed in Sec. IIVI where it is shown to be in agreement with Edwards' 
thermodynamic description. The compactivity is identified in terms of the parameters char- 
acterizing the mesoscopic dynamics of the model. Also, the distribution of domains is there 
derived. Section contains a detailed discussion of the relationship between closed and 
open models, and between the compactivity and fugacity parameters. The paper ends with 
a short summary and some additional discussions. 

II. THE MODEL 

We consider a one- dimensional lattice having N + 1 particles. The number of sites in 
the lattice is variable and changes in time accordingly with the rules to be specified below. 
Those sites that are not occupied by a particle are said to be empty or, equivalently, being 
a hole. 

The dynamics of the system is defined trying to mimic tapping experiments for the 
study of compaction in granular media These experiments typically involve two 

different series of processes of quite different nature. The system is submitted to taps or 
pulses separated by time intervals for which the system is allowed to relax freely, until 
being trapped in a metastable configuration. Therefore, each tap starts in the metastable 
configuration reached in the previous free relaxation. The taps are characterized by their 
duration and their amplitude. 

Physically, the effect of the taps is to decrease the local density in some regions of the 
system, moving grains from their metastable positions, and allowing a posterior reordering in 
the free relaxation. We will specify first the dynamics during the relaxation processes, since 
it leads to identify the possible metastable, or blocked, configurations of the model. It will 
be assumed that in the free relaxation the system tries to reduce its length by eliminating 
some of the empty sites of the lattice. More precisely, whenever there is a group of nearest 
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neighbour holes, all except one are eliminated. These are the only processes taking place in 
the free relaxation, and have probability one. Therefore, the number of particles is conserved 
in the relaxation, but the length of the system, measured by the total number of sites, is in 
general reduced. 

As a consequence of the above dynamics for relaxation, the metastable configurations of 
the model are characterized by having all the holes surrounded by two particles, i.e. the 
holes are isolated. In order to displace the system from one of these states, it has to be 
externally perturbed, for instance by means of a tap. To complete the description of the 
dynamics of our model in a compaction experiment, the possible transitions taking place 
during a tap and starting at a mestastable configuration, have to be identified and their 
probabilities specified. Two kinds of elementary processes will be considered. Each of them 
will be discussed separately in the following. 

Firstly, a particle can be transiently desorbed from the lattice and posteriorly adsorbed 
in an empty site in the neighbourhood of its previous position. This process is restricted 
by the following rule. A particle can be desorbed from a site during the tap only if at least 
one of its nearest neighbour sites is empty. More precisely, the probability for these events 
is proportional to the number of nearest neighbour holes of the particle being desorbed. 
This restriction tries to naively model the short range constraints making difficult structural 
rearrangements in granular materials. Then, during a tap, the probability of desorption of 
a particle having only one nearest neighbour hole is a, while it is 2a if it is surrounded by 
two holes. Afterwards, the particle is reabsorbed either in its own original site or in any of 
the nearest neighbour holes, with a probability that is proportional to the number of holes 
next to the site considered. In Fig. ^ cases (a) and (b) involve processes starting with the 
transitory desorption of a particle. Particles and holes are represented by circles and crosses, 
respectively. In the case referred to as (b) in the figure, the elimination of a hole happening 
in the next free evolution has been also indicated. It is seen that the net result of the series 
of events taking place during the tap and the free relaxation is, in this case, the destruction 
of a lattice site, with the consequent decrease of the lattice length. 

During the tap, the creation of an empty site or hole is also possible, but only between two 
nearest neighbour particles located at one of the ends of a domain of at least two particles. 
The probability of the corresponding elementary events, referred to as case (c) in Fig. Q 
is (3. Note that these processes of hole creation are just the inverse of those producing the 
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FIG. 1: Elementary rearrangements of the system in a single tap and the following free relaxation, 
in the weak vibration limit. The trajectories leading to a final state identical to the initial one are 
not shown. 

destruction of a hole. 

Furthermore, it will be assumed that only one transition takes place at the most in every 
region on the system during each tap, i.e. no site is involved into two different processes 
in the same perturbation of the system. Physically, this hypothesis implies to consider the 
limit of weak and short taps 0,0]. In summary, the dynamics of the model in the shaking 
process is defined by the effective transitions given in Table HI describing the combined events 
associated to a tap and the next free relaxation. The transitions only affect the clusters 
shown, and their probabilities are independent of the configuration of the remainder of the 
system. 

To formulate the model in a more mathematical language, and also to characterize the 
metastable configurations, it is convenient to define a set of variables n = {m, • • - , njv}- 
The variable takes the value unity if there is a hole next to the right of particle i, while 
it vanishes if there is no hole, i.e. if particles i and i + 1 are in nearest neighbour sites. By 
definition, it is assumed that there is no hole to the left of particle 1 nor to the right of 
particle N + 1. Both particles define the boundaries of the system. It is easily realized that 
this property is preserved by the dynamics of the system under tapping, as defined above. 
Then, we have established a one to one relationship between a set of N variables taking 
values and 1 and the metastable configurations of the model. 

The transition probabilities in Table |l] can be expressed in terms of the variables. 
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TABLE I: Transition probabilities for the elementary rearrangements taking place in a single tap, 
in the weak vibration regime. 



Initial state Final state Probability 

OOXO 0X00 a/2 

0X00 00X0 a/2 

0X0X0 0X00 a/2 

0X0X0 00X0 a/2 

0X00 0X0X0 p 

00X0 0X0X0 p 



Denoting i?jn = {ni, • • • , RiUi, ■ ■ ■ , njsr}, with i^rij = 1 — nj, the probability W(n'|n) of 
the several events going from configuration n to configuration n' in the effective dynamics 
describing shaking process are: 

Oi 

W(RiR i+1 n\n) = - [m(l - n i+1 ) + (1 - ni)n i+1 ] , (1) 
oi 

W{Ria\a) = -(rii-ini + nin i+1 ) + f3 [ni_i(l - ni) + (1 - ni)n i+1 ] . (2) 



Equation (JTJ corresponds to diffusive events, while the first and second terms on the rhs of 
Eq. (J2J) correspond to the destruction and creation of holes, respectively. The Markov process 
defined by these transition probabilities is irreducible, i.e. all the metastable configurations 
of the lattice are connected by a chain of transitions with nonzero probability. To verify 
this property, we begin by noting that any metastable configuration can be connected with 
the configuration characterized by having just a hole located next to the right of a given 
fixed particle. This is because holes can be moved through particles by means of diffusive 
events, so that two consecutive holes can always be located to both sides of the same particle. 
Afterwards, one of the holes can be eliminated by a type (b) event of Fig. This procedure 
can be repeated until there is only a hole in the lattice, that can then be diffused to the 
desired site. This proves the above statement. But, since each effective transition have its 
inverse also with nonzero probability, the above paths can also be reversed, concluding that 
all the metastable configurations are connected. The irreducibility property of the Markov 
process implies that there is a unique steady probability distribution for the process 
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This distribution will be explicitly obtained in the next Section. 

In Fig. El the relaxation of the particle density is shown as a function of the "scaled time" 
t = an, where n is the number of taps before measuring the density, for different values of 
a and (3. The initial state for all the curves was the least dense metastable configuration, 
p = 0.5, in which there is a hole between every two particles. In all the reported cases 

C a, so that processes decreasing the density of particles are only relevant when the 
system is near the most compact state, p — 1. Moreover, as (3 <C a <C 1, there is a universal 
behaviour up to a very large number of taps t = 0{a~ 1 ) 1 i.e., an = 0(1). For longer times, 
when processes decreasing the density become relevant, n > 0(f3~ l ), the system approaches 
a steady state characterized by the ratio (3 /a. The observed universal scaled curve is very 
well described by the four-parameter empirical law 

pit) = poo c- , (3) 

l+51n 1 + ^ 



with poo = 1.04, 5poo = 0.54, B = 1.17 and r c = 2.63. As it is the case with the experimental 
data pj and also with numerical results from other simple models jiol lisj]. the logarithmic 
fit is not expected to give the correct asymptotic density of particles. In fact, in our case it 
is Poo > 1, which is clearly unphysical. A similar behaviour of poo was found in Ref. jlol |. 
Also, values of poo larger than the random close packing limit have been reported from the 
fit of experimental data j]]. 



III. THE STEADY DISTRIBUTION 



To find the steady distribution of the Markov process describing the effective dynamics 
of the model, we are going to suppose it verifies detailed balance. Of course, this has to be 
justified a posteriori by showing that such a distribution exists. Therefore, we look for a 
time-independent distribution p( s \n) having the property 

iy(n'|n)p (s) (n) = W(n|n> (s) (n') , (4) 

for all configurations n and n'. A direct first consequence of this equation is that all the 
metastable configurations with the same number of holes have the same probability in the 
steady state. This follows from the fact that they are connected by diffusive events and 
diffusion is isotropic in the effective dynamics, as seen in Table H] Therefore, the distribution 
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FIG. 2: Evolution of the density of particles, as a function of the scaled time defined in the text. 
The curves correspond to the pairs of values a = 10" 3 , (3 = 1(T 5 (circles), a = 10" 2 , (3 = 5 x KT 5 
(squares) and a = 10 -2 , j3 = 5 x (diamonds), while the solid line is the best inverse logarithmic 
fit, Eq. PJ ). with the parameters given in the text. 



function verifying Eq. Q can only depend on the number of holes, 



N 



ri; 



in the configuration, but not on their spatial distribution. So, we can write 



(5) 



(6) 



where the number iV# of holes in the configuration n has been made explicit in the notation, 
f(Nn) is a function to be determined, and Z a normalization constant, 

Z = J2f(N H ). (7) 

n 

Still remains to be analyzed if Eq. (jUJ) can be made compatible with Eq. (J3J), when par- 
ticularized for effective events increasing (and decreasing) the number of holes. The latter 
reads 

(8) 



Here n' is a configuration differing from n by the creation of a hole. Use of Eq. © gives 

f(N H + 1) _ 2(3 
f(N H ) a 



(9) 
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and, by iteration, 

/(•V//) rl-l , (io) 



for Njj > 1, C being an arbitrary constant that will be taken equal to unity. In this way, 
we have proven that the system has the property of detailed balance and that its steady 
distribution is given by 

= (11) 

N 

z = Y,t Nh = E n N {N H )>y- N *, (12) 

n N H =1 

where 7 = a/2/3 and £In{Nh) is the number of metastable configurations of the lattice having 
Nh holes and, of course, N + 1 particles. It is worth remarking that no approximation has 
been done in order to derive the steady distribution, Eq. (JTTJ, i.e, it is valid for any value 
of 7. The steady average number of holes and its dispersion can be evaluated from Z by 

( Ns ). = Z"*Ato = ( 13 ) 

n J 

A simple combinatorial argument gives 

= w (J £W (15) 

and substitution of this expression into Eq. (|12|) yields 

N 

7. 

Therefore, it follows by application of Eq. (|13p that, in the limit of large N 



1 + 1) -1. (16) 



W. = y^. (17) 

The right hand side of Eq. (fT7|) is a monotonic decreasing function of 7, for fixed number of 
particles N, i.e., the length of the system decreases as 7 increases. Therefore, 7 -1 plays in the 
model a role similar to the vibration intensity in real granular experiments of compaction. It 
follows that, in the physical image depicted by the present model, the probability of diffusion 
proceses a is expected to grow faster with the vibration intensity than the probability of 
creation of holes /3, al least in the weak tapping limit. 
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FIG. 3: Comparison between the numerical values of the steady density of particles and the 
theoretical prediction, given by Eq. lfTU|) . 



The length (volume) of a configuration is given by 



L = N H + N 



(18) 



and the average particle density is 



P 



N 



N 



1 



(19) 



(L) s (N H ) S + N 2 + 7 

In Fig. El the numerical values for the steady density of particles, obtained by Monte Carlo 
simulation of the model, are compared with the theoretical prediction, Eq. (|19jl. and an 
excellent agreements is found. The specific length per particle, in site units, is the inverse 
of the particle density, 



<0« 



N+(N H ) S 2 + 7 



N 1+7 
Its dispersion is obtained from Eqs. ()14j) and ()16p. 

iV((A/) 2 ) s = (2-(/) s )((/) s -l), 



(20) 



(21) 



presenting a maximum for (l) s = 3/2, i.e. when the average number of holes is N/2 and the 
density of particles = 2/3, i.e., 7 = 1. The numerical evaluation of the length fluctuations 
is compared with the theoretical prediction, as given by Eq. (|21j). in Fig. 0] Again, a very 
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FIG. 4: Comparison between the numerical evaluation of the length fluctuations and the theoretical 
prediction, given by Eq. (|!?Tjl. 

good agreement is observed for the range of "vibration intensities" 7 -1 plotted. Outside this 
window of vibration intensities, the length fluctuations are very small and, therefore, rather 
hard to measure in the simulations. 



IV. THERMODYNAMIC DESCRIPTION 



Following Edwards and coworkers ideas |4J , the steady distribution (fTTj) can be expressed 
in the "canonical" form 

-N H /X 

= Z = ^n(N h )c- n ^ x , (22) 

where 

X = (1U7)" 1 (23) 

is the so-called compactivity. It is the conjugated thermodynamic parameter of the volume in 
gently vibrated granular systems, in an analogous way as the temperature is the conjugate of 
the energy in usual thermal systems. Note that, in Eq. f!22[) . the ratio Nh/X can be replaced 
in both the numerator and the denominator by L/X, where L is the length (volume) of the 
configuration, as defined in Eq. The structure of the above steady distribution is 



12 



consistent with the two main ingredients of Edwards' theory, namely that the measure 
over metastable configurations is flat, and that there is a unique parameter, the volume, 
characterizing the macroscopic state of the system. Let us point out that, very recently, the 
theory has been extended to include several macroscopic control parameters, in an effort 



12, lia, and 



to explain the discrepancies observed in some models with strong tapping 
also segregation patterns in binary models Q]. It is clear that such an extension does 
not apply to our model, that is designed to describe compaction in one-component systems 
under weak tapping. In the same context, the expression of the compactivity in Eq. ()23j) 
deserves some comments. Although X can be formally negative, for values 7 < 1, it is quite 
doubtful that this fact be physically relevant, since this range of values of 7 corresponds to 
strong tapping, leading to low stationary densities, namely with a average number of holes 
(Nh) s > N/2. The possibility of negative values of the compactivity has been also found in 
other simple models jloL 20], and it is associated with the existence of a maximum length 



for the metastable configurations. 

In the limit 7 — > 00, i.e., asymptotically weak tapping, the steady concentration of 
particles, p( s \ given by Eq. (JIHJ), can be approximated by 

p (s) ~ 1 - - (24) 
7 

and, using the definition in Eq. (J23J), 

X- 1 = -ln(l -p (s) ). (25) 

This relation between the compactivity and the steady density has been also found in a 
model with facilitated dynamics having a variable number of particles and fixed volume 
[ic| . and a similar behaviour has been reported from the analysis of experimental data jsj]. 
An "entropy" S associated to the distribution can be defined in the usual way, 

S = -5>W(n)lnp (s) (n), (26) 

n 

and use of Eq. (J2*2~j) gives 

s = iM 1 + lnZ (27) 

X 



Taking into account Eq. (|13jh it is easily verified that 

as as 1 

diT)~s ~ d(N H ) s ~ X' 
13 



(28) 



consistently with the physical meaning of the compactivity as discussed above. Given that 
the macroscopic state of the system is characterized by a single parameter, it is possible to 
express the entropy in terms of only the density of particles, or the intensity parameter 7, 
or the compactivity. Then, for instance, in the limit of large N the entropy can be written 

as 

S _ 1 1 + e 1 ^ 

N~ X{l + e^ x ) + e V* • ( 9) 

In addition to the global properties considered up to now, it is also possible to obtain 

information about the domain structure of the steady configurations. In particular, we are 

going to derive here the probability distribution for the number of particles in a domain. 

A domain of size r is defined as a cluster of r particles, i.e. two holes with r particles in 

between. First, we consider the probability of finding a local domain of size r, 

Fj: s) = (n fc (l - nfc+i) • • • (1 - n fc+r _i)n fc+r ) s , (30) 

with r < 1. Use of Eqs. and (|15|) yields 

1 N-r-l f \ r+1 

F { r s) = ^ E ^WA^h" 2 -^ = 7" 2 M- > ( 31 ) 
z Nh=o V 1 + 7/ 

where the limit of large N has been considered once again. Then, the probability of a domain 
of size r, P(r), is given by the conditional probability of finding a cluster of r consecutive 
particles plus a hole to the right of a given hole, i.e., 

p v = wriw&=v+w- (32) 

This distribution is correctly normalized: 

00 

£P(r) = l. (33) 

r=l 

It is instructive to express the distribution of domain sizes in terms of the average length 
per particle, (l) s . This is easily accomplished by means of Eq. (|2())1. obtaining 

P(r) = (2 - (O.r 1 ((0, - 1) ■ (34) 



V. RELATIONSHIP BETWEEN CLOSED AND OPEN MODELS FOR COM- 
PACTION 

In the previous Section, we have introduced the compactivity X from the canonical form 
of the steady probability distribution, Eq. (|2*2*|) . In the Edwards and coworkers formulation 

14 



of the granular thermodynamic theory |4j, the compactivity was defined by 

'dS\ 

dv) 



N 



where the entropy S is given by 

S = \nn N , (36) 

Qn being the number of blocked configurations or, in the language used in this paper, 
metastable states. In Eq. (|35|) . the number of particles in the system is kept constant. The 
quantity Qn for the model considered in this paper is given by Eq. (|15p. and for N ^> 1, 
N H > 1, Eq. flUSJ) leads to 

*-> = l»^. (37) 

This is the "microcanonical" (constant volume) version of the "canonical" (constant com- 
pactivity) expressions (|17j) and (|23J1 . In fact, combination of these two latter expressions 
gives 

x - 1 = ln A^. (3g) 

The expression equivalent to Eq. ()36|) in the canonical ensemble is Eq. f!26|) . It is evident 
that, in the limit of large systems, it is consistent with the definition of X in Eq. ([35)1 . 

In several proposed models for compaction, lattices with a fixed number of sites, i.e. 
fixed length, have been considered. The dynamics is defined involving elementary processes 
associated with the adsorption and desorption of particles, in such a way that the number 
of particles in the lattice changes along the shaking experiment. This is the mechanism 
for which the density in the system varies in time. In particular, several models leading to 
similar kind of metasta 
in detail 
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e configurations as in the model in this paper have been discussed 



14L Il9| . Then, aside from details that are irrelevant for the following 



analysis, the number of blocked configurations is given by Eq. (J15)) . that we rewrite in the 
form 

= N^L-W H )V < 39 > 
where the number of sites of the lattice L is now considered as fixed and L/2 > Nh > 1. 
Moreover the steady distribution, in the weak tapping limit was found to have the form 

P (s) \n) = V —, (40) 
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with 

L/2 

Z> = £ n N {N B )rT N *, (41) 

N H =1 

and r) is a given parameter, depending on the specific model, and characterizing the dynam- 
ical events in the system under shake. Then, from Eq. (J4*U|) a compactivity X' was identified 

as 

X' = {lnr])-\ (42) 

This definition is equivalent to 

dn^iv^ fas' 

or, using Eq. (JHHJ), 



^ _ ^ ~ , (43) 



This expression differs from Eq. (J37)l except in the limit of high density iV# /./V <C 1, in which 
both reduce to \ti(N/Nh), But this only indicates that the same density is obtained in this 
limit if X = X', although it must be stressed that X is associated to tapping processes at 
constant number of particles, while X' describes processes at constant volume. Equivalently, 
X characterizes ensembles with fixed N, and X' ensembles with fixed L. In this context, 
their physical nature is rather different. The parameter X is the compactivity introduced 
by Edwards and, on the other hand, related with X' by Eq. (|42|) . plays the role of a 
fugacity for the particles. In terms of the entropy, X and X' are related by 

X- l = X*-*+(%L) • (45) 



ON 



N H 



VI. CONCLUSION 

In this paper, a one-dimensional model for compaction in granular media has been pre- 
sented. One of its main features, as compared with previous Ising-like models, is that the 
time evolution under tapping conserves the number of particles in the system, while it is the 
volume what changes in the compaction process. This is in fact what happens in compaction 
experiments. Consequently, the steady distribution is characterized by the compactivity in- 
stead of by a generalized fugacity. The steady distribution function has been derived and 
the compactivity identified in terms of the parameters defining the mesoscopic dynamics of 
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the model. It has been found that the results are consistent with Edwards' thermodynam- 
ical theory of powders. Nevertheless, since the model is formulated in the context of weak 
and short tapping, it is in fact quite doubtful that the same conclusion were reached from 
a generalization to stronger tapping processes. Let us point out that this would require to 
modify the formulation of our model by including the possibility that a lattice region would 
experiment several elementary excitations during the same tap. 

The relationship between closed and open models, and between compactivity and fugacity, 
has been discussed. At the mesoscopic level of description used in this paper, the expression 
of one of them in terms of the transition rates can not be inferred from the expression of 
the other. Nevertheless, it is true that they correspond to different derivatives of the same 
entropy function, like in usual thermal systems. 

The model presented here can be easily generalized to mixtures of several kinds of grains, 
then allowing the study of segregation phenomena. Also, it can be useful to investigate the 
validity of the Edwards theory in this case, and eventually its possible generalizations, for 
instance by extending the number of parameters needed to characterize the steady state of 
the mixture, as it has been proposed recently [lfl llfi|. 
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